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Defects in the honeycomb lattice of graphene are known to induce local moments and strong 
correlation effects. Distortions due to structural reconstruction around vacancies in graphene were 
studied recently by Cazalilla et al. [arXiv: 1207.3135 (2012)], who formulated an effective model 
consisting of a localized a level hybridized with the 7r-band. We analyze the rich quantum impurity 
physics of this system using a combination of numerical renormalization group and analytical tech- 
niques, focusing on the special role played by the unusual local density of states, which is enhanced 
at low energies due to potential scattering. Depending on microscopic parameters, the model hosts 
both exactly-screened spin-1/2 (doublet) Kondo or underscreened spin-1 (triplet) Kondo phases, 
and we study the quantum phase transition separating them. Although the effective Kondo models 
also support new stable phases characterized by strong renormalized particle-hole asymmetry, such 
phases cannot in fact be accessed in the full Andersonian model describing the vacancy. We show 
that distinctive signatures of the modified powerlaw Kondo effect thus always appear at low energies 
in thermodynamic quantities and the scattering t matrix. 

PACS numbers: 73.22.Pr, 72.15.Qm, 72.10.Fk 



Since its experimental isolation, graphene has at- 
tracted a great deal of attention on both experimental 
and theoretical fronts alike. 1 3 Remarkably, most of the 
novel properties of graphene can be explained within a 
simple non- interacting single particle picture. However, 
richer physics naturally results when strong electron cor- 
relations and formation of magnetic moments play a key 
role, as is the case for example when transition metal 
atoms are adsorbed onto the graphene surface. 

In standard metals, the Kondo effect describes screen- 
ing of such local moments by quantum many-body spin- 
singlet formation at low temperatures, 4 which constitutes 
a classic paradigm of strong correlations in condensed 
matter science. But due to the peculiar nature of the 
low-energy degrees of freedom in bulk graphene, 5 10 in- 
cluding the vanishing density of states at the Dirac point, 
magnetic impurities may exhibit unusual Kondo physics, 
and even local quantum phase transitions. 10 14 

Defects in graphene, such as carbon vacancies and 
induced reconstructions in the otherwise perfect hon- 
eycomb lattice, can also give rise to interesting new 
physics. 15 19 In particular, evidence of the Kondo effect 
has been observed experimentally in irradiated graphene 
samples which host such vacancies. 20 22 It was argued 
recently by Cazalilla et al. in Ref. 19 that structural cor- 
rugations around vacancies allow for a hybridization be- 
tween the a- and 7r-orbitals, which leads to local magnetic 
moment formation (an effect absent in flat graphene). 
Potential scattering from the defect induces a greatly en- 
hanced density of conduction electron states coupling to 
these local moments. 19 This can potentially give rise to 
very high Kondo temperatures, as observed already in 
certain experiments. 22 For a thorough recent review of 
Kondo physics in graphene, see Ref. 23. 

In this paper we study in detail the effective model for 



isolated vacancies in graphene given in Ref. 19. The very 
unusual local density of states leads to a variant of the 
pseudogap Kondo problem, giving rise to exotic quantum 
impurity physics. For example, the residual impurity en- 
tropy of Si mp = — &b In (4) is the lowest possible for this 
type of system. The complex orbital structure of the 
underlying impurity problem also supports a quantum 
phase transition separating doublet and underscreened 
triplet Kondo phases. We go beyond the previous pertur- 
bative analysis, using the full numerical renormalization 
group (NRG) to calculate various physical quantities ex- 
actly (for a recent review of the technique, see Ref. 24). 
Numerical results are supported wherever possible by an- 
alytical arguments. 

Our main results are as follows: 

(i) Phase diagrams for the effective Kondo models with 
spin S = 1/2 (doublet) and S = 1 (triplet) are presented 
and discussed. In both cases, the modified powerlaw 
density of states supports two distinct strong coupling 
phases, corresponding to runaway RG flow of either the 
Kondo exchange coupling or potential scattering. Ther- 
modynamic and dynamical quantities are analyzed ana- 
lytically at the various stable fixed points, while the full 
crossover between fixed points is calculated numerically. 
The quantum critical physics associated with the phase 
transition separating strong coupling states is studied for 
both doublet and triplet Kondo models. 

(ii) The quantum phase transition between S = 1/2 and 
S = 1 Kondo screened phases of the full Andersonian 
model is analyzed. We find that the quantum phase tran- 
sition is a simple level crossing. 

(iii) We find that the particle-hole asymmetric local mo- 
ment phases of the simplified Kondo models are not ac- 
cessible in the underlying full Andersonian model, and 
hence will not in practice be realized in experiments. 
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The organization of the manuscript is as follows: Fol- 
lowing Ref. 19, we introduce the coupled two-impurity 
Anderson model describing the reconstructed graphene 
vacancy in Sec. I. Depending on the strength of the 
Hund's rule coupling, both S =1/2 and S = 1 'impurity' 
configurations can be realized. Deep in either regime, ef- 
fective low-energy Kondo models can be derived, which 
are studied explicitly in Sec. II. Thermodynamic and dy- 
namical quantities are calculated in each phase of each 
model, and their phase diagrams are analyzed in detail. 
Physical behavior in the vicinity of the stable strong cou- 
pling fixed points is understood and interpreted in terms 
of a resonant level or pure potential scatterer. The evo- 
lution of strong coupling scales is extracted in each case. 
We also examine the quantum phase transition separat- 
ing Kondo screened and asymmetric local moment states 
in the doublet and triplet Kondo models. 

In Sec. Ill we analyze the full Andersonian model, and 
argue that certain phases of the Kondo model are not in 
fact accessible. However, a level-crossing quantum phase 
transition separating doublet and triplet Kondo screened 
phases is accessible, and we study it in detail. We con- 
clude with a discussion of our results in Sec. IV. More 
technical parts of the analysis are relegated to the ap- 
pendix. 



I. MODEL AND OBSERVABLES 

The underlying Anderson model for a single vacancy 
in graphene was presented in Ref. 19. It is formulated 
in terms of two hybridized local quantum orbitals of the 
a and 7r bands, denoted d a and 7r a respectively. The 
spin density of these orbitals, = \S a a (7( jid (j , and 
= \t\) t (j ( j ( j<t\ (jI , interact via the Hund's rule cou- 
pling Jh- Local electron correlations on the d-level, 
parametrized by J7, and a capacitive interaction U an act- 
ing between d and tt levels are also important. 19 We al- 
low for tunneling of electrons between the levels as well 
as the hybridization of the d-level with the conduction 
electrons of the graphene host. The full Hamiltonian is 
H = H + #im P + H hyh , where 

k,cr 

#im P = e d n d + Unln^ + e n n n + U an n d n n 

-JtfS d -S^+#^(47T a +H.C.) , 

a 

H hyh = v hyh J2 ( c k,<A + H.c.) . (1) 

k,cr 

Here, c^ a denotes the conduction electrons, and n d = 
5^cr n d = 5^cr d^dp is the total number operator for the 
d-level, while for the 7r-level, = n° = ^\^ a ^ 

Importantly, the local conduction electron density at 
the defect site, p(uo) oc ^ k S(uj — ek), is strongly modified 



by potential scattering. 19 It assumes the diverging but 
normalizable form 

^ \uj\ m \oj/A\ 

where Do is a high-energy cutoff for the conduction band 
(the precise choice of which is irrelevant for the low- 
energy quantum impurity physics). We stress that the 
second cutoff A > Do is necessary for the density of states 
to be normalizable. In a real system, A corresponds phys- 
ically to the scale above which the density of states is not 
well-described by Eq. (2) because high-energy details of 
the graphene band structure become important. In par- 
ticular, we note that the Kondo temperature Tk will thus 
naturally depend on the couplings and also on A. 

The retarded free local electronic Green function at 
the defect site follows from the local density of states, 
Eq. (2), as 

G iocM = V [de^-Z7rp(u), 
J (jJ — e 

^ ln ^ (2 oj \ 
2|w|ln 2 |^| \ir A J 

(3) 

where V denotes the principal value, and the second line 
is obtained analytically at low frequencies, as sketched in 
Appendix A. 

The hybridization function for the d- level, A (a;) = 

^hyb^ioc^)' con tains all the information about the non- 
interacting bath of conduction electrons and, together 
with i^i m p, completely specifies the underlying quantum 
impurity problem. As demonstrated in the following, 
the model can then be solved numerically exactly within 
NRG by exploiting a logarithmic discretization of this 
hybridization function. 

For realistic model parameters, 25 there are two rele- 
vant regimes of the isolated d-it level system described 
by i^imp- The parameter G = essentially controls 

the occupancy of the 7r-level. For G <C 1, the 7r-level 
is essentially empty, (n^) ~ 0; while for G ^> 1, the n- 
level is singly-occupied, (n n ) ~ 1. When the d-level is 
also singly-occupied (as obtained for — U <C e d <C 0), 
the ground state of i^imp is thus a doublet for G <C 1, 
but a triplet when G ^> 1 by virtue of the ferromagnetic 
Hund's coupling. 



A. Effective low-energy Kondo models 

On hybridization with the graphene conduction elec- 
trons, low-energy effective Kondo models can be derived 
using the standard Schrieffer- Wolff transformation, 4 pro- 
vided v^ yh /U <C 1. Projecting onto the doublet or triplet 
manifolds of i^imp by perturbatively eliminating virtual 
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charge fluctuations to second order in Uh y b> one obtains 
H K = H + J K S • s + V J2 c k, CT c k , )CT (4) 

k,k',cr 



where S 
and So : 



: Sd + is the total spin of the d-ir subsystem; 



k,k' 



7 /c k , a ,. In the doublet regime 
we thus have an S = 1/2 Kondo model, while in the 
triplet regime an S = 1 Kondo model results. Since 
Kondo models with different impurity spin give different 
ground states, 4 the low-energy mapping of the full model 
Eq. (1) onto the effective Kondo models Eq. (4) imply 
the existence of a quantum phase transition on tuning 
G. This scenario is explored in detail in Sec. IIIB. The 
effective parameters of these Kondo models are given by, 



J K ~ 2|v hyb | 5 

y „ Kbl 2 



1 



1 



1 



1 



\U + e d \ 



(5) 



where e d = e d in the doublet case, while in the triplet case 
£d = e d + U a7T — Jh/^- It is important to note that in both 
cases the Kondo coupling Jk > is antiferromagnetic. 
Of course, the novel feature of these Kondo models is 
the unusual density of states, Eq. (2), whose modified 
powerlaw behavior is shown to result in rich quantum 
impurity physics. 



B. Physical observables 

In the following, we consider thermodynamic proper- 
ties of the above Andersonian and Kondo models. In 
particular, we wish to calculate impurity contributions 
to the entropy, Si mp (T), and the magnetic susceptibility 

Ximp 

(T) = ((£tot) 2 )im P /2\ as a function of temperature, 
T. As usual, 4 (f2)i mp = (Q) — (O)o? with (Cl)o denoting 
the thermal average for the free system. The behavior 
of these quantities allows straightforward identification 
of the various fixed points, and evince clear RG flow. 
Crossovers between these fixed points also provides di- 
rect access to the underlying energy scales of the prob- 
lem, such as the Kondo temperature, Tk- 

In addition, we study dynamic quantities, focusing 
on the energy-dependence of the t matrix. Since the t 
matrix controls the scanning tunneling spectroscopic re- 
sponse and resistivity measurements, it is the key quan- 
tity needed to interpret certain experimental results. The 
rapid and universal evolution of the t matrix at energies 
on the order of the Kondo temperature could thus pro- 
vide the 'smoking gun' signature of Kondo physics in a 
graphene system with isolated vacancies. 

The t matrix describes scattering between diagonal 
states k and k' of Ho induced by the defect. It is gener- 
ically given by, 



(6) 



where 7k,k'(^) are components of the t matrix. Since 
hybridization between the d-level and the host is local in 
space, one readily obtains, 



G loc (") = < ) M+ GEM T{u), 



: (0), 



(7) 



where T(uo) = NTk^(uj) and N — » oo is the number of 
graphene conduction orbitals. 

For the full Anderson model, Eq. (1) the t matrix fol- 
lows as 



7k,k'M = vl yb G dd (u), 



(8) 



where G dd {uo) = ((d a ; * s ^ ne d-\eve\ Green function. 

At low energies in the Kondo regime, the same scat- 
tering should be produced by the corresponding Kondo 
models, Eq. (4). The t matrix is then expressed as, 



where 



V 



(9) 



(10a) 



Jk/2 



((7;7 f ))c , (10b) 



with 7 = S z coa + S C04 (here c , a = ^ Ek c k,<x is a 
local operator for the conduction electrons at the defect 
site and S z , S~ denote operators for the impurity spin). 
7jf k' i 00 ) 18 tne trivial contribution to scattering from a lo- 
cal potential, while 7^\/ (w) describes the effect on scat- 
tering due to electron correlations, such as the Kondo 
effect. Below, we consider explicitly the spectrum of the 
t matrix, i.e., 



t(u) = Im T(u) 

7T 



-N Im 7L, k ' • 



(11) 



II. KONDO MODELS 



As discussed above, when charge fluctuations are 
frozen out on the d- and ir- levels, the low-energy physics 
of the reconstructed vacancy is described by a Kondo 
model, with either an effective S = 1/2 or S = 1 local 
moment coupled to the graphene conduction electrons. 
These models are a close relative of the pseudogap Kondo 
problem, 11 14 ' 26 30 where the density of states is given 
generically by p(uj) oc |co>| r . Such models have been stud- 
ied extensively, especially for S = 1/2 in the context of 
certain high-T c superconductors and magnetic impurities 
in regular graphene, 10 which both realize the r = 1 case. 
Vojta and Bulla considered the case — 1 < r < 0, de- 
scribing a spin- 1/2 impurity coupled to a diverging den- 
sity of states. 31 Although their survey allowed to identify 
the topology of the phase diagram, details of the various 
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crossovers and properties of the fixed points themselves 
were not established exactly. Further detailed study of 
the generalized powerlaw Kondo model for arbitrary an- 
tiferromagnetic or ferromagnetic Kondo coupling will be 
considered in a separate publication. 32 



A. Analytical results 

1. Generalized poor man's scaling 

Before solving the model numerically exactly using 
NRG, we first consider perturbative scaling for the S = 
1/2 and S = 1 Kondo models, Eq. (4), with density of 
states given by Eq. (2). A similar analysis was performed 
in Ref. 19 for a Kondo model with a related density of 
states, although the role of the cutoff A was not appreci- 
ated. 

We demonstrate below that two scaling invariants arise 
for both S = 1/2 and S = 1 models, corresponding to 
RG flow of the Kondo coupling Jk and potential scat- 
tering V. Strong coupling scales Tk and Tp can thus be 
estimated, and are shown to depend on A as well as the 
bare couplings and the band cutoff Do. 

In the spirit of Anderson's poor man's scaling 19 ' 31 ' 33 
we derive flow equations in terms of the dimensionless 
couplings j = D £Z R In ^ and v = D ^ 2 » In ^, viz: 



dj 



d\nD 

dv 
d\nD 



= -1- 



(12) 



Remarkably, the second equation describing scaling of 
the potential scattering term is exact to all orders. The 
coefficient of the linear terms is the scaling dimension 
and here it depends upon the running scale D and the 
bare 'cutoff' A. Both scaling equations can be solved 
analytically to yield the scale dependent couplings: 



3(D) 
v(D) 



P{D) 



3(D ) 



p(D ) i _i_ MPs. 
P(D) 



p(D a ) Sd P( D' )dD ' 



P(D ) 



v(Do) 



(13) 



An immediate observation is that the strong coupling 
phases vanish as A — » Dq since then j(D), v(D) —> 0. 
The density of states, Eq. (2), reduces to a pair of poles 
at ±D in this limit, and so there are no low-energy bath 
degrees of freedom to which the impurity spin can cou- 
ple. However, in the generic case A > Do, the running 
Kondo coupling j(D) and potential scattering v(D) are 
renormalized upward on successive reduction of the band- 
width/energy scale. The strong coupling scales Tk and 
Tp associated with the divergence can be identified ap- 
proximately as the point where these couplings become 



of order unity: i.e. j (Tk) = 1 and v (Tp) = 1. Of 
course, such a perturbative approach breaks down before 
this point, but such an analysis does provide a useful 
analytic estimate. 

To leading order, it suffices to consider the scaling di- 
mension of the flow equations. To this level of approx- 
imation, the flow of j and v is identical, and so the re- 
sulting strong coupling scales have the same functional 
dependence, 



T P In 



2 

A 

2 



Jk, A 
— In , 

2 D ' 

V. A 
2 ln IV 



(14) 



In the limit of small Jk (V), and consequently small Tk 
(Tp), Eq. (14) can be inverted to yield 



T K 



Jk lnj 
2 ln 2 ^ 



2 In' 



v 

A 



(15) 



The notable feature here is that the strong coupling scales 
are proportional to the bare couplings Jk and V. Of 
course, this is in marked contrast to the standard metal- 
lic case, where p(uj) = po is finite and constant at low 
energies. An exponentially small Kondo temperature, 4 
Tk oc exp(— 1/pqJk), results in that case. Remarkably, 
strong coupling physics in the defective graphene system 
should be expected at relatively high temperatures. 19 

Below, we extract the exact dependence of Tp on V in 
the pure potential scattering case (Jk = 0). In Sec. II C 
we go further and extract Tk and Tp from exact non- 
perturbative NRG calculations for the non-trivial case 
where V ^ and Jk ^ 0. 

The implication from the RG scaling equations, 
Eq. (13), is that at high energies the physics is controlled 
by a 'local moment' (LM) fixed point, describing a free 
impurity spin, decoupled from the bare host (the impu- 
rity spin is either S =1/2 or 5=1, as above). For 
larger bare Jk (^> V), the Kondo coupling is renormal- 
ized faster as the temperature/energy scale is reduced, 
and hence the Kondo effect will dominate. For 5 = 1/2, 
the ground state is then described by a 'symmetric strong 
coupling' (SSC) fixed point, and the impurity spin is 
screened by conduction electrons. The remaining Fermi 
liquid host degrees of freedom now also feel an additional 
phase shift (equivalent to a modified boundary condition 
at the defect site) due to the Kondo effect. In the present 
case of a diverging free density of states, this results in 
an additional anomalous host contribution to physical 
properties, as shown below. For impurity spin S = 1 the 
impurity is in fact under screened, 4,19 and the fixed point 
is denoted 'USSC'. Alternatively, if the potential scatter- 
ing is initially very strong (V ^> Jk), then V is renor- 
malized faster and the ground state is described by the 
'asymmetric local moment' (ALM) fixed point. Although 
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the impurity (either spin S =1/2 or 5 = 1) remains free, 
the renormalized potential scattering in the host becomes 
maximal and similarly generates an anomalous contribu- 
tion to physical properties. A critical point, (U)AF-CR, 
separates the (U)SSC and ALM phases. The topology of 
the flow diagram is shown in Fig. 1, and explicitly dis- 
cussed in Sec. II B where the physical picture is confirmed 
by means of exact numerics. The underlying topology of 
the phase diagram is thus equivalent to that studied in 
the pure powerlaw case by Vojta and Bulla. 31 

2. Potential scatterer 



ujp thus serves as a definition of the strong coupling scale, 
Tp = up. The exact result in the limit Jk = is then, 

T P ln(T P /A) = In A , ( 20 ) 

which should be contrasted with the perturbative scaling 
result, Eq. (14). The two definitions will be compared in 
Sec. (II C). 

The change in thermodynamic quantities due to the 
introduction of such potential scattering can also be ob- 
tained from the electron Green function, Eq. (6), once 
the t matrix is known. 4 At low temperatures, the behav- 
ior in the vicinity of the ALM fixed point is found to be, 



The simplest limit of the Kondo model arises for 
Jk = in Eq. (4). The model is then trivial in the sense 
that it is non-interacting, and Kondo physics per se is 
totally absent. However, as suggested by the flow equa- 
tion, Eq. (13), this simple system has a strong coupling 
ground state, characterized by potential scattering whose 
strength diverges below an emergent scale, Tp. Later, it 
will be shown that the low-energy behavior of the Kondo 
model with Jk > can also be understood in terms of 
the pure potential scatterer throughout the ALM phases. 
We discuss briefly this limit now. 

The local Green function at the defect site is modified 
due to the additional potential scattering for V ^ 0. It 
is expressed exactly as, 

-l 



[GiocMr 



V 



(16) 



Rearranging this equation in the form of Eq. (7) immedi- 
ately yields the t matrix for the pure potential scatterer, 
Eq. (10a). The corresponding spectrum has the following 
asymptotic forms, 



t(u) = 



( 1 i i i I s S n M i i 1 2 i 00 

+ OQoj\ 3 ), (17a) 
1 



Do V|w|ln 2 |w/A|, 

+ 0(1/M 2 ). (17b) 

By demanding dt(u)/du = 0, we find that the spectral 
peak for small V arises when, 



Vp = —up In 



Up 



/In 



A_ 
A) ' 



(18) 



As up — >> 0, the spectrum exhibits scaling in terms of 
u/up, viz: 



t(u) |w|<|w P | 



u 



Up 

|2- 




sgn[co>V] 

+ 0{\u/u P \ 2 ) , (19a) 
0(\w p /uj\ 3 ) . (19b) 



S imp (T) = -ln(2) 1 + 



Inf. 



TXimp(T) = 



Inf. 



+ 0(T 2 ) , (21a) 



+ 0{T 2 ) . (21b) 



As confirmed by full NRG calculations below, this leading 
behavior arises on the lowest energy scales in the ALM 
phase also in the case of finite Kondo coupling, Jk > 0. 



3. Resonant level 

In the Kondo screened SSC phase of the pseudogap 
Kondo model, it has long been established 13 ' 14 ' 27 that 
the lowest-energy physics is that of a simple resonant 
level model, 



H RL =H +vY,(cl a d a +H.c.) 



(22) 



k,cr 



which is just the non- interacting (U = 0) Anderson im- 
purity model at particle-hole symmetry. Indeed, in the 
present case of the defect in graphene, we show below 
that the lowest-energy behavior of the SSC phases can 
again be understood in terms of the resonant level. This 
motivates a brief analysis of the latter for the case where 
H is given by Eq. (1), with corresponding local density 
of states, Eq. (2). 

The retarded resonant level Green function is simply, 



1 



u-v 2 G{° o l(uj)+i0+ 



(23) 



where G[^ c (u) is given in Eq. (3). The corresponding 
spectrum is given asymptotically by, 



t{u) 



M«H / 1 



1 



2 In 



v 2 In 2 A. 



I 3 In 



+ 0(M 5 ), (24a) 



l~4i A ( 1 

= kv In — — k 

2 Do VM 3 m 2 |w/A| 



+ 0(1/M 5 ). (24b) 



The position of the spectral peak at ujk is straightfor- 
wardly obtained, and the spectrum again obeys scaling 
as uok in terms of uj/ujk. 



t(uj) \oj\^[uj k \ 



1 



1 



w|»|wjf| 








U 2 t) 





+ 0(\uj/uj k \ a ) , (25a) 
+ 0(KH 6 ), (25b) 



where uj 2 K hi\uj K /K\ = -£ 2 ln(A/D ). 

Hrl should be regarded as an effective low-energy 
model here, valid in the vicinity of the SSC fixed 
point, 13 ' 14 and so the parameter v is itself an effective 
parameter. Thus, one cannot directly identify ujk for the 
resonant level model with Tk for the full Kondo model 
with Jk > 0. However, the asymptotic scaling form of 
the spectrum, Eq. (25a), is expected to hold at low ener- 
gies, as confirmed in the next section. 

To lowest order, the RG scaling equations given in 
Eq. (13) are the same for the pure potential scatterer 
(with J K = 0) and the Kondo model (with V = 0). The 
strong coupling scale for the former is given by Eq. (20); 
mutadis mutandis, the leading dependence of the Kondo 
temperature in the latter should be, 



T K ln(T K /A) = -\J K \\n 



Do 



(26) 



On the lowest energy scales in the vicinity of the SSC 
fixed point, thermodynamics of the full spin- 1/2 Kondo 
model, Eq. (4), can be calculated from the effective res- 
onant level model. 13 ' 14 ' 27 The leading low-temperature 
behavior can again be extracted exactly from the elec- 
tron Green function (or the corresponding resonant level 
t matrix), as shown explicitly for the entropy in Appendix 
B. As T -> we find, 



Simp(T) = -lll(4) 1 + 



1 



+ (9(T 2 ) , (27a) 



TXimp(T) 



1 1 + — U I +(9(T 2 ) . (27b) 



In 



2T 
A , 



The low-temperature form of Eq. (27) is confirmed below 
by explicit NRG calculations for the full S = 1/2 Kondo 
model in the SSC phase. 



B. Numerical results 

The physics of the spin-^ Kondo model with finite Jk 
and V [Eq. (4)] is obviously much more rich and subtle 
than the trivial limits considered above. Here, one gener- 
ically expects two phases: for Jk/\V\ ^> 1 a strong cou- 
pling (U)SSC phase should result; while an ALM phase is 
stable for \V\/Jk ^> 1- These phases are separated by a 




LM 



oo 
(U)SSC 



FIG. 1: Schematic phase diagram for the Kondo model with 
density of states given by Eq. (2) and 'impurity' spin S =1/2 
or S = 1. \V\ denotes particle-hole symmetry breaking while 
Jk is the Kondo coupling. Fixed points are denoted by circles, 
while arrowed lines refer to effective RG flow. The location 
of the separatrix (dashed line) depends on spin S, although 
the topology of the diagram is unaltered. 



quantum critical point ((U)AF-CR) arising for a critical 
ratio (J K /V) C = 0(1). 

The full temperature-dependence of thermodynamic 
quantities can be calculated using the NRG technique. 24 
Their characteristic behavior at the various fixed points 
allows straightforward identification of the phases, and 
the entire phase diagram can thus be mapped out. We 
find that the topology of the phase diagram for spin 
S = 1/2 and antiferromagnetic Jk > is the same 
as for the case of the pure powerlaw density of states 
(p(u) ~ \uj\ r with -1 < r < 0) studied in Ref. 31. In- 
deed, the topology is unchanged for S = 1, although the 
fixed points themselves describe different physics and give 
rise to different thermodynamic signatures. A schematic 
phase diagram for the Kondo model with density of states 
given by Eq. (2) and spin 5 = 1/2 or S = 1 is presented 
in Fig. 1 and discussed further below. 



1. S = 1/2 case 

At high energies/temperatures, the LM fixed point de- 
scribes a free and unscreened impurity local moment. 
The limiting high-temperature entropy is thus Si mp = 
ln(2), while the magnetic susceptibility follows the Curie 
law, Tximp = 1/4. But the Kondo effect drives the sys- 
tem toward the SSC fixed point below an energy scale 
of order T K when J K /\V\ > (J K /\V\) C . For T <C Tk 
the entropy and susceptibility for the spin- 1/2 Kondo 
model are found from NRG to be given precisely by 
Eq. (27), obtained for the effective resonant level model. 
In particular, we note that the residual T = entropy 
Simp(0) = — ln(4) is the smallest possible value, because 



7 



S = l/2 S = l 




T/D T/D 

FIG. 2: Thermodynamics across the quantum phase tran- 
sition for S = 1/2 (left panels) and S = 1 (right pan- 
els). Entropy Si mp (T) (upper panels) and magnetic suscep- 
tibility TximpiT) (lower panels) plotted vs T/Do for fixed 
A/Do = 1-5, Jk = 0.1A) and tuning V — ► K- Shown for 
V = F c =b 10- n T%=° with n = 0,1,2,3,4,5 in order of de- 
creasing T* (oc \V — Vc\). Solid lines for the (U)SSC phase; 
dashed lines for the ALM phase. The critical point itself is 
shown as the dotted lines. For comparison: diamond points 
for Jk/Dq =0.1, V/Do = and circle points for Jk/Dq = 0, 
V/D = 0.1. 



the conduction electron density of states is characterized 
at low energies by the strongest possible divergence (up 
to logarithms), while remaining normalizable. 

By contrast, for J K /\V\ < (Jk/\V\) c the ALM fixed 
point is stable, and describes maximal particle-hole 
asymmetry. The crossover from LM to ALM physics oc- 
curs on the strong coupling scale, Tp. The free impurity 
local moment is then supplemented by an anomalous host 
contribution to give T — » thermodynamics described 
by Eq. (27). The precise agreement on the lowest energy 
scales confirms the physical interpretation of the fixed 
point in terms of the pure potential scatterer discussed 
above. 

On fine-tuning in the vicinity of the critical point 
Jk/\V\ ~ {Jk/\V\) c , RG flow from the high-energy LM 
fixed point to the stable fixed point describing the ground 
state occurs via the critical point, AF-CR. Two universal 
scales can thus be identified in this regime: T c (« Tk) 
sets the scale for onset of criticality associated with AF- 
CR, while T* (oc \V — V c \) characterizes the ultimate 
crossover to either SSC or ALM fixed points (depend- 
ing on the sign of V — V c ). The distinctive asymptotic 
thermodynamic properties of each fixed point are sum- 
marized in Table I. 

Remarkably, the impurity entropy at ALM, AF-CR, 
and SSC fixed points is negative. We note, however, 
that all thermodynamic quantities flow to their respective 



TABLE I: Properties of fixed points for 5=1/2 





lmiT^O 'S'imp 


lim T ^o Tximp 


ALM 


-In 2 


1/8 


LM 


In 2 


1/4 


AF-CR 


-In 4/3 


1/24 


SSC 


-In 4 


-1/8 



zero-temperature values logarithmically slowly. This is a 
direct consequence of the logarithmic energy dependence 
in the density of states, see Eq. (2) 

The left panels of Fig. 2 show the full crossovers in 
the entropy and susceptibility calculated numerically ex- 
actly using NRG. In both cases, diamond points show the 
direct crossover from LM to SSC arising for Jk/Dq = 
0.1 but V/Dq = 0; while circle points show the di- 
rect crossover from LM to ALM when V/Dq = 0.1 but 
Jk/Dq = 0. Thermodynamics are also shown on tuning 
to the quantum critical point (dotted line). Approach 
from the SSC phase (solid lines) and the ALM phase 
(dashed lines) exhibit two characteristic scales, T c and 
T*, as above. 

Dynamical quantities, such as the T = scattering 
t matrix, similarly evince the rich RG structure of the 
problem. In Fig. 3 (upper panel) we plot the scaling 
spectra t(u)/\u\ vs \uo/uok\ for a system deep in the SSC 
phase (solid and dotted lines). Its asymptotic behavior 
at both high and low energies is found to be described 
by Eq. (25). 

In the ALM phase (Fig. 3, center panel), the scaling 
spectrum t(uu)/\uj\ vs \u/up\ is asymptotically described 
by Eq. (19). In particular, we note the more gentle linear 
approach to the Fermi level, and the inherent asymmetry 
of the spectrum, arising due to the relevance of particle- 
hole symmetry-breaking. 

The lower panel of Fig. 3 shows a critical spectrum, 
plotted as t(uj)/\uj\ vs \uj/uj c \. The resonance around uo = 
sgn[cj]/]|cj c | is split, with the peaks separated by ~ \V — 
Jk\- At low energies |cj| <C \uo c \ the spectrum t(uo)/\uo\ ~ 
a + b sgn[ujV]\uj/uj c \ has a leading linear dependence. 



2. S = 1 case 

In the case of a spin-1 Kondo 'impurity' (obtained at 
low-energies from Eq. (1) when G = ^rf- ^ 1)? the 
phase diagram is topologically equivalent to that of the 
spin S = 1/2 case (see Fig. 1). The thermodynamic 
properties of the fixed points themselves are summarized 
in Table. II. 

TABLE II: Properties of fixed points for S = 1 





limT^O Simp 


lim T ^o Tximp 


ALM 


-In 4/3 


13/24 


LM 


In 3 


2/3 


UAF-CR 


In 5/4 


3/8 


USSC 


-In 2 


1/8 



8 




\w/ UK,P,c I 

FIG. 3: Spectral function t(uj)/\uj\ at T — 0. For impurity 
spin S = 1/2: cj > and uj < plotted as solid and dotted 
lines. For impurity spin S = 1: cj > and cj < plotted 
as dashed and dot-dashed lines. Top panel: vs \uj/ujk\ in 
the SSC phase (J K = 10 _6 A) and V/Jk = and 0.01). 
Middle panel: vs \uj/ujp\ in the ALM phase (V = 10~ 6 D and 
Jk/V = and 0.01). Lower panel: vs \uj/uj c \ at the critical 
point (J K = 10~ 6 D and V/J K = (V/J K )c). A/A) = 1.5 
used throughout. Asymptotes discussed in the text. 



The LM fixed point again describes the high- 
energy /temperature behavior of the 5=1 Kondo model. 
However, the free impurity moment is now a spin-1 ob- 
ject, which gives rise to 5i mp = ln(3) and Tximp = |, as 
seen clearly for all data in the right-hand panels of Fig. 2. 

For Jr/\V\ < (Jk/\V\)ci relevant particle-hole sym- 
metry breaking drives the system to the ALM fixed point 
below Tp. The free spin-1 impurity moment is then sup- 
plemented by the anomalous host contribution to give 
S'imp = m (3) — ln(4) and Tximp = § — §■ The direct 
crossover from LM to ALM is shown as the circle points 
in the right-hand panels of Fig. 2. 

For Jk/\V\ > (Jk/\V\) c , the Kondo effect drives the 
system to the USSC fixed point below Tk- For spin-1, 
the impurity is however under screened: the conduction 
electrons only partially compensate the impurity mo- 
ment. Asymptotically, a free impurity spin- 1/2 moment 



remains. Together with the anomalous host contribution, 
one then obtains 5i mp = ln(2) — ln(4) and Tximp = \ — §■ 
The direct crossover from LM to USSC is shown as the 
diamond points in the right-hand panels of Fig. 2. 

In the vicinity of the quantum critical point, Jk/\V\ — 
(Jk/\V\) c , two scales again emerge. Flow toward the un- 
stable UAF-CR critical point on the scale of T c (~ Tk) 
is followed ultimately by flow to the ALM or USSC fixed 
points below T* (~ \V — V c \), which describe the T = 
ground state of either phase. The right-hand panels 
of Fig. 2 also show the full temperature-dependence of 
the impurity contribution to the entropy, 5i mp (T) (up- 
per panel), and magnetic susceptibility, Tximp(^) (lower 
panel), as the critical point (dotted line) is approached 
from the USSC phase (solid lines) and the ALM phase 
(dashed lines). As with the spin- 1/2 case, the approach 
to all fixed points occurs logarithmically slowly as a di- 
rect consequence of the explicit logarithmic energy de- 
pendence of the density of states, Eq. (2). 

Scaling spectra t(uo)/\uo\ at T = for systems in the 
USSC phase (upper panel, vs \u)/wk\), ALM phase (mid- 
dle panel, vs \u/up\) and at the UAF-CR critical point 
(lower panel, vs \uo/uo c \) are shown in Fig. 3 as the dashed 
and dotted lines. They are characterized by the same 
asymptotic behavior as the spin-1/2 case, described by 
Eqs. (19) and (25). 



C. Evolution of T K and T P 

Solution of the generalized poor man's scaling equa- 
tions for the spin 5=1/2 and S = 1 models gives an es- 
timate for the strong coupling scales Tk and Tp, Eq. (14). 
In the trivial potential scattering limit considered explic- 
itly in Sec. II A 2, the scale Tp can be obtained exactly, 
Eq. (20). Up to logarithms, one thus expects in either 
case the same linear dependence of the strong coupling 
scale Tp ~ V or Tk ~ Jk- 

Here we calculate Tp and Tk numerically exactly using 
NRG 34 for a wide range of systems with 5 = 1/2 and 
5 = 1, varying V, Jk and the cutoff A. In Fig. 4 we 
plot Tk ln(A/7V)/ ln(A/£)o) vs Jk/D$ for systems in the 
(U)SSC phase, and T P ln(A/T P )/ ln(A/D ) vs V/D for 
systems in the ALM phase. The excellent agreement, 
especially for small Jk or V, confirms Eqs. (26) and (20) 
(solid line). For comparison, the result of Eq. (14) is 
shown as the dashed line. 



III. ANDERSON MODEL 

We consider now the full Andersonian model for the 
graphene vacancy, Eq. (1). 
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10" 6 10" 4 10" 2 10° 

Jk/Do, V/D 

FIG. 4: Tk and Tp extracted from NRG calculations for var- 
ious values of Jk and V deep in the SSC phase (Jk /V — 2) 
and the ALM phase (Jk/V = 1/2) for both 5 = 1/2 and 
S = 1 and A/ D = 100, 10, 1.1, 1.01, 1.001. All data plotted 
as TK,pln(A/T K ,p)a| ) p/ln(A/ J Do) vs J K /D or V/D [the 
constant a KP = 0(1) depends only on the phase (SSC or 
ALM) and on the impurity spin S] . Data collapse to the ana- 
lytic result [Eq. (20), solid line], is found at small Jk, V. For 
comparison, Eq. (14) for A = 10 is shown as the dashed line. 




Jk/D 

FIG. 5: Phase diagram in the Jk — V plane for A/ Do = 1.5, 
indicating the critical separatrix between (U)SSC and ALM 
phases for impurity spin S =1/2 and S = 1. At small Jk, the 
critical ratio \V/Jk\c is given by Eq. (29) [see inset, dashed 
lines]. The regime of validity of the effective Kondo model 
is given by the Schrieffer- Wolff asymptote \V/ Jxlmax = 1/4, 
shown as the dotted line. 



A. Accessibility of the ALM phases 

In Sec. II we discussed the rich phase diagram of the 
Kondo models with 5 = 1/2 and 5 = 1. In both cases, 
these models were derived from an Andersonian model 
(which also allows for charge fluctuations), by project- 
ing onto the 'impurity' spin-manifold using a Schrieffer- 
Wolff transformation. Within this leading-order pertur- 
bative treatment, Eq. (5) indicates that the maximum 
ratio of the effective potential scattering, V, and the ef- 
fective Kondo exchange coupling, Jk is given by 

V 1 

* £ 5- (28) 

The natural question is then: which of the strong cou- 
pling phases of a Kondo model can actually be accessed 
within the bare Anderson model. In order to answer this, 
we study the exact position of the phase boundary of the 
effective 5=1/2 and 5 = 1 Kondo models; and examine 
directly the full Andersonian model, Eq. (1). 

Fig. 5 shows the phase diagram of the Kondo models 
for both spin 5 = 1/2 (diamond points) and 5 = 1 (circle 
points), obtained numerically exactly using NRG. The 
inset shows the asymptotic small- Jk behavior, which is 
found to follow 

=a 5 + 6In(J^/A), (29) 

with a s =\/2 — 3/4, as=i = 1 and b ~ —0.3 for A/D = 
1.5 (see dashed lines). 

For comparison, we also plot the Schrieffer- Wolff re- 
sult, Eq. (28), as the dotted line in Fig. 5. The obvious 



conclusion is that the effective potential scattering, V, 
derived from the bare Andersonian model is never large 
enough to access the ALM phase. A numerical survey of 
the parameter space of Eq. (1) supports this result, and 
suggests that the Kondo effect is always operative at the 
lowest energy scales. The ground state is thus described 
by the SSC (or USSC) fixed point. 

We note, however, that this conclusion strictly applies 
only to the model with simplified density of states given 
in Eq. (2). Very strong particle-hole asymmetry in the 
density of states at high energies could in principle change 
this picture. 

B. The doublet-triplet transition 

Tuning the Hund's coupling, J#, in Eq. (1) allows to 
access either 5 = 1/2 or 5 = 1 sectors of the isolated 
'impurity' manifold. On coupling to the conduction elec- 
trons for finite Vh y t>, effective low-energy Kondo models 
can then be derived with definite impurity spin, which 
have different stable fixed points (SSC for 5 = 1/2 but 
USSC for 5 = 1). Since the two ground states cannot 
be continuously connected (they have, for example, dif- 
ferent residual entropy), one expects a quantum phase 
transition on tuning Jh in the full model, which we call 
DT-CR. 

This scenario is explored in Fig. 6, where we take rep- 
resentative model parameters for the Andersonian model, 
and fine-tune Jh across the transition. Solid lines cor- 
respond to the SSC doublet phase [Jh < (Jh)c], while 
dashed lines are for systems in the USSC triplet phase 
[Jh > (Jh)c\- The critical point itself is shown as the 
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FIG. 6: Thermodynamics 5 f imp (T) and T X im P (T) vs T/D 
across the quantum phase transition, separating doublet and 
triplet phases of the full model, Eq. (1), on tuning the 
Hund's coupling Jh- Parameters are A/ Do = 1.5, U/Do = 
0.3, e/D = -0.15, U a7T /D = 0.025, v hyh /D = 0.008, 
g/D = 0.008, varying J H /D = (J H )c ± 10~ n T^ H=0 with 
n = 0, 1, 2, 3, 4 approaching progressively the critical point at 
(J H )c ^ 0.101A) (dotted line). For J H < (Jh)c the doublet 
phase is obtained (solid lines); while Jh > (Jh)c yields the 
triplet phase (dashed lines). The transition is characteristic 
of a level-crossing. 



TABLE III: Properties of fixed points for the full model 





limT-)-0 Smp 


lhTlT^O ^Ximp 


ssc 


-ln4 


-1/8 


ussc 


-ln2 


1/8 


DT-CR 


-In 4/3 


1/24 



dotted line. The asymptotic thermodynamic properties 
of the fixed points are also given in Table III. 

We find that the doublet-triplet transition is charac- 
teristic of a level-crossing: the low-energy scale setting 
the crossover from DT-CR to either SSC or USSC has 
the linear dependence T* ~ | Jh — (Jh)c\- Physical prop- 
erties at the point of transition can be understood simply 
in terms of coexisting SSC and USSC fixed points. 

Deep in either doublet SSC or triplet USSC phase, we 
have confirmed by direct numerical calculation that the 
low-energy physics is that of the corresponding 5 = 1/2 
or S = 1 Kondo models in the SSC or USSC phases. 



IV. CONCLUSIONS 

In this paper we studied the physics of a reconstructed 
vacancy in an otherwise perfect graphene sheet. The de- 
fect is modelled in terms of a generalized Anderson model 
with a diverging density of states following a recent paper 
by Cazalilla et al. 19 We combined analytical and numeri- 
cal techniques to obtain a rather complete understanding 
of the model and its properties. 

When charge fluctuations are frozen out, the com- 
plex orbital structure of the Andersonian model can be 
more simply understood in terms of an effective spin- 1/2 
or spin-1 object exchange-coupled to the graphene host. 
These doublet and triplet Kondo models are unusual due 
to the modified graphene conduction electron density of 
states, which results in a strongly enhanced linear scal- 
ing of the Kondo temperature with coupling strength, 
Tk ~ Jk (up to logarithms) and anomalous ground state 
thermodynamic properties. For example, in the doublet 
phase, the residual impurity contribution to the entropy 
assumes the minimum possible value for a quantum im- 
purity problem of this type, S{ mp (T —*►()) = —feeing). 
Furthermore, it leads to a rich phase diagram, because 
particle-hole symmetry-breaking is relevant. Unlike the 
situation in standard metals, a new stable strong cou- 
pling phase is supported in these Kondo models, charac- 
terized by large renormalized potential scattering, which 
suppresses the Kondo effect. A more conventional strong 
coupling phase is also supported, associated with a large 
renormalized Kondo coupling (although anomalous host 
contributions to certain physical quantities are again 
shown to appear). The two phases are separated by an 
unstable quantum critical point. We supplement our de- 
tailed numerical study with analytic treatment of simple 
effective models which describe the physics in the vicinity 
of the two strong coupling fixed points. 

However, direct analysis of the full Andersonian model 
reveals that the phase transition separating the two 
strong coupling phases cannot in practice be accessed, 
because the effective Kondo coupling and effective po- 
tential scattering are slaved. In consequence, the Kondo 
effect is always operative on the lowest energy scales in 
the graphene system with isolated vacancies, modelled 
by Eq. (1). In principle, potential scattering from other 
sources could manifest at higher energies (\uj\ > A, not 
considered here). If sufficiently strong, we note that the 
particle-hole asymmetric local moment phase could be 
stabilized; although the distortion of the density of states 
required for this makes it unlikely. 

The existence of two effective Kondo models with 
distinct spin implies the presence of a quantum phase 
transition in the full Andersonian model between 
exactly-screened spin- 1/2 Kondo and underscreened 
spin-1 Kondo phases. We demonstrated that this 
transition can be realized on tuning the Hund's cou- 
pling, and find that the transition is characteristic of 
a first-order level crossing. We study both phases and 
the transition between them using the NRG technique. 
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As expected from the SchriefTer- Wolff mapping, the 
low-energy physics of the full model for the isolated 
graphene vacancy is indeed shown to be described by 
effective Kondo models in their Kondo screened phases. 
The unusual physics of the generalized Kondo problem 
with strongly divergent conduction density of states is 
thus expected in such systems. 
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This expression directly leads to Eq. (3). 



Appendix B: Logarithmic corrections 



Appendix A: Local Green function 

The derivation of the local Green function at low fre- 
quencies is straightforward but tedious. In the following 
we sketch it: while the imaginary part follows trivially 
from the density of states the real part has to be evalu- 
ated from Kramers-Kronig relation 
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uj — e 
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The real part can be brought into the form 
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After an integration by parts one arrives at 
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At this point it is important to notice that in the limit 
\uj\ <C A we can rewrite the above expression to leading 



Here we obtain analytically the first logarithmic cor- 
rection to the impurity entropy, from the free energy of 
an effective resonant level model. The general formula to 
calculate the free energy of a local level is given by 



(Bl) 



Using residual calculus this can be converted into a line 
integral along the branch cut along the real axis given by 

/CO J 
^n F (z)e^ 0+ Imln(-G-i(z)) (B2) 
-oo 7r 

where the subscript r refers to the retarded Green func- 
tion. Consequently, we obtain the entropy as 



S= - 
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(B3) 



where the factor of two is due to the spin summation. 
The low-energy properties will be entirely dominated by 
the self-energy, which is of the form 
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It is then straightforward to derive an expression for the 
impurity entropy. It follows as, 
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Realizing the low-temperature limit, this can be approx- 
imated as, 
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which, to leading order, yields 
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